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(Dated: August 5, 2008) 

We present a study of B decays into semileptonic final states containing charged and neutral 
Di(2420) and D2(2460). The analysis is based on a data sample of 208 fb"^ collected at the 7(45) 
resonance with the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. With a 
simultaneous fit to four different decay chains, the semileptonic branching fractions are extracted 
from measurements of the mass difference Am — m{D**) —m{D) distributions. Product branching 
fractions are determined to be B{B+ ^ D^^i^e) x S(L>? D*+tv~) = (2.97 ± 0.17 ± 0.i7) x 10~^ 
B{B+ Dfi+ue) X B{Df D'^*'>+tt~) = (2.29±0.23±0.21) x 10-^ B{B° D-£+vi) x B{D- 
D*°n-) = (2.78±0.24±0.25) x 10"^ and B{B° ^ D;-£+iye) x S(D;" D'-*^°-k-) = (1.77±0.26± 
0.11) x 10"^. In addition we measure the branching ratio r(Z)2 Dn )/r(Z)2 D^*'tx ) = 
0.62 ±0.03 ±0.02. 

PACS numbers: 13.20.He, 13.25.Ft, 14.40.Lb 



Measurements of the Cabbibo-Kobayashi-Maskawa 
matrix elements \ Vcb\ and \ Vub\ rely on precise knowledge 
of semileptonic i?-meson decays. Decays with orbitally 
excited charm mesons (£)**) in the final state give a sig- 
nificant contribution to the total semileptonic decay rate. 
A better understanding of these decays will reduce the 
uncertainty in the composition of the signal and back- 
grounds for inclusive and exclusive measurements [1]. 

In the framework of Heavy Quark Symmetry (HQS), 
D** mesons form two doublets with — 1/2^ and 
= 3/2~ where denotes the spin-parity of the light 
quark coupled to the orbital angular momentum. The 
doublet with = 3/2~, namely the Di and D^, have to 
decay via D wave to conserve parity and angular momen- 
tum and therefore are narrow with widths of order of 10 
MeV 01 ■ The relative contribution of the two doublets 
and the polarization of the produced D** mesons can be 
compared with QCD sum rules |3] and predictions from 
Heavy Quark Effective Theory [4j . 

In this Letter we describe a simultaneous measurement 
of all B semileptonic decays to the two narrow orbitally- 
excited charmed states, without explicit reconstruction of 
the rest of the event. The CLEO collaboration has pre- 
viously reported a branching fraction measurement for 
B+ Dli+v and an upper limit for B+ Df£+u 5]. 
Belle and BaBAR^ have reported results using a technique 
in which one of the B mesons in the process T{4S) BB 
is fully reconstructed 0. 

In this analysis we use a sample with a total integrated 
luminosity of 208 fb~^, part of the complete data set col- 
lected with the BaBAR detector at the PEP-II storage 
ring, operating at a center of mass energy of 10.58 GeV. 

The BaBAR detector Q and event reconstruction fsl are 
described in detail elsewhere. A Monte Carlo (MC) sim- 
ulation of the detector based on GEANT4 is used to 
estimate signal efficiencies and to understand the back- 
grounds. The sample of simulated BB events is equiv- 
alent to approximately three times the data sample and 
a dedicated simulation of signal events based on the 
ISGW2 model |10] has been produced with statistics 
equivalent to roughly 5 times the expected signal yield 
contained in the data. 



D** decays are reconstructed in the decay chains 
D**^D*TT- [m, and D**^Dtt~. The former is acces- 
sible to both narrow D** states while the latter has no 
contribution from the Di. Intermediate D* states are 
reconstructed in D*—yD'^Tr and the D mesons are recon- 
structed exclusively in D'^^K~Tr~^ and D'^^K'tt'^tt'^ . 
D** candidates are then paired with reconstructed lep- 
tons and required to be consistent with the semileptonic 
decays B — > D**£v, as described in the following. 

First, events which are most likely to contain a semilep- 
tonic B decay are selected. We require that there is a 
reconstructed D candidate and at least one lepton in the 
event with a momentum greater than 800 MeV/c [l^. D° 
meson candidates are formed by K~ti^ combinations re- 
quiring the invariant mass to be consistent with the 
mass: 1.846 < m{KTi) < 1.877 GeV/c^. This asymmet- 
ric mass window is chosen to take into account resolution 
effects of the detector. The selection is optimized to max- 
imize the significance of the selected sample. 

D'^ candidates are combined with charged and neu- 
tral pions to form D* candidates. For D*^ the 7r° is 
reconstructed from a photon pair with an invariant mass 
of 115 < TO-y^ < 150MeV/c^. Those photon pairs are 
re-fitted in a "mass-constrained" fit to match the nomi- 
nal mass of the 7r°. D* candidates are selected by their 
mass difference to the candidate: 144 < m{D^TT^) — 
m{D°) < 148MeV/c2 and 140 < m{D°TT") - m{D°) < 
144 MeV/ c'^ for charged and neutral D* , respectively. 

candidates are formed from K^tt'^tt^ combina- 
tions with an invariant mass of 1.854 < m{K'K'K) < 
1.884 GeV/c^ The fit probability for the three tracks 
to originate from a common vertex, Pytx, is required to 
be Pyt^iKinr) > 0.01. 

Candidates for D and D* are combined with charged 
pions to form D** candidates, and finally paired with 
muons or electrons. The charge of the lepton is required 
to match the charge of the kaon from the D decay. 

Part of the background is due to events where a D** is 
paired to a lepton from the other B. Thus we require that 
the probability that the lepton and the pion emitted by 
the D** originate from a common vertex exceeds 0.001, 
and that the angle between the direction of flight of the 
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D** and the lepton is more than 90 degrees. 

A large fraction of the background events is due to 
B D*lv decays where the D* or its daughter D is 
paired to a pion from the other B. To suppress this 
combinatorial background, we make use of the variable 
COSBY described in the following. The energy and mo- 
mentum of the B mesons from the T{AS) decays are 
known from incident beam energies. For correctly re- 
constructed B D**(.v decays, where the only miss- 
ing particle is the neutrino, the decay kinematics can be 
calculated, up to one angular quantity, from the four- 
momentum of the visible decay products {Y ~ D**£). 
The cosine of the angle between the direction of flight of 
the B meson and its visible decay product Y is given by 

2EbEy — m^B ^ 



where E, \p\ and m are the energies, momenta and masses 
of the B and the Y, respectively. If the Y candidate is 
not from a correctly reconstructed B ~> D**lv decay, the 
quantity cos by no longer represents an angle, and can 
take any value. We select candidates having jcossy | < 1. 

In case a D* is reconstructed in the decay chain, a 
veto is applied against decays B D*£i' by calculating 
the variable cos by ' which is defined as above, but the 
Y system is redefined to contain only the D* and the 
lepton: Y' = D*i. Background events are rejected by the 
requirement cos by' < — 1 since signal events B — > D**tv 
tend to have values less than —1. 

To reduce combinatorial backgrounds in the decay 
chain D**^Dtt^ , only the D**l candidate with ml 
closest to zero is selected, where is the neutrino 
mass squared, calculated in the approximation ps = 0: 
™J = "^B + l-Pi'P ~ 2EbEy- Events reconstructed in 
the D**^D'^tt^ final state are rejected if the 1?° can be 
paired with any charged pion to form a candidate 
as described above. 

In about 2% of the events more than one D**£ candi- 
date is selected and if so all of them enter the analysis. 

We determine the signal yield in the channel 
D** Dtt and the Di and D2 signal yields in the 
channel D** — > D*^ by a binned fit to the Am = 
m{D^*\') — ■m{D^) distributions. To determine the in- 
dividual contributions from Di and D2 in the D*tt final 
state, we make use of the helicity angle distribution of the 
D*, -dh, which is defined as the angle between the two pi- 
ous emitted by the D** and the D* in the rest frame of 
the D* . For a D* from a D2 this distribution varies as 
sin^ "dh, whereas for Di decays, the helicity angle is dis- 
tributed like 1 + cos^ -dh, where Ab^ is a parameter 
which depends on the initial polarization of the Di and a 
possible S wave contribution to the Di decay. To exploit 
this feature, we split the data for the two decay chains 
involving a D* into four subsamples, corresponding to 
four equal size bins in | cos "dh \ ■ 



The resulting ten Am distributions are fitted simulta- 
neously to determine 12 parameters describing the sig- 
nal yields and distributions, and 22 parameters to ad- 
just the background yields and shapes. The mass dif- 
ferences for the signal events are described by Breit- 
Wigner functions. There are four parameters giving the 
signal yields for the semileptonic decays involving the 
two narrow states, charged and neutral. The masses 
of the states are also fitted, but are constrained to be 
equal for charged and neutral states, giving two pa- 
rameters. Four additional parameters arise from the 
effective widths of the D** states, which represent a 
convolution of the intrinsic widths and detector reso- 
lution effects. The latter contributes approximately 2- 
3MeV/c^, depending on the mode. The fit also de- 
termines the D2 branching ratio Bb/dm = ^{D2 
Dtt-)/ {r{D* ^ Dtt-) + r{D* ^ D*t:-)) and the Di 
polarization amplitude Abi- 

Backgrounds are modeled by cubic functions in Am. 
The background shape in the D*Tr~ channel is found to 
be the same in all helicity bins for each final state. The 
fit thus has three shape parameters for each decay chain, 
while the number of background events is determined 
independently in each bin. 

The selection efficiency is deduced from a fit to the 
simulation. This fit uses the same parametrization as 
the fit determining the signal yield from data and is ap- 
plied to the sum of the full background simulation and 
for one signal decay chain at a time. For a given de- 
cay mode the efficiencies are found to be the same for 
Di and D^, specifically: e(i:>*+7r-) = (6.89 ± 0.12)%, 
e{D*°7T-) = (5.34±0.12)%, e{D+7r-) ^ (12.88 ± 0.96)% 
and e{D°Tr-) = (17.56 ± 0.70)%, where the quoted un- 
certainties are the statistical uncertainties from the fit. 
For the decays including a D* the efficiency is mul- 
tiplied by the probability for a D** to decay with a 
value of IcosiJ/il falling into a given bin. This fac- 
tor includes the theoretical distribution discussed above 
as well as corrections for the different detector accep- 
tances in the four helicity bins of up to 10%. The total 
number of B mesons in the data sample used for the 
present work is Ngg = (236.0 ± 2.6) x 10^ [13]. For the 
charged and neutral B mesons we assume TiTjAS) 
B+B-)/r(r(45) B°B°) = 1.065 ± 0.026 Q. 

The fit procedure has been extensively validated. The 
analysis procedure is tested on statistically independent 
MC simulated data samples and was found to reproduce 
the input signal parameters with a x^/?i = 12.66/12, 
where n is the number of signal parameters. Consistent 
fit results were also obtained when the data sample was 
separated into subsamples representing specific data tak- 
ing periods, separated by lepton species or restricting it 
to certain decay modes, using charged or neutral D** 
only, or combining the helicity bins. The results of the 
fit are shown in Fig. [TJ As expected, the contribution 
of the D2 vanishes for large values of | cos \ while the 
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Am (GeV/c-) Am (GeV/c^) 



FIG. 1; Am-spectra for the selected data and the results of the fitted functions. The solid line represents the complete fit 
function, dotted (Di) and dashed (DJ) lines for the signal and dash-dotted the for background, (a) to (d) show the mode 
D**°^D'+iT- with increasing values for | cos (e) the mode D**'^—>D'^n . (f) to (i) show the corresponding bins in | cos '&h\ 
for the mode D"+ -*D*'\+ and (k) the mode D"+ ^D°tt+ . 



contribution of the Di is suppressed for cos'dh close to 
zero. The extracted yields are given in Table [H 

Systematic uncertainties have been analyzed and their 
impact on the fitted yields have been estimated taking 
into account correlations between fit parameters. Effi- 
ciencies for reconstructing and selecting the particles of 
the final state are derived from Monte Carlo simulation. 
The simulation of the tracking and the tt" reconstruction 
have been studied by comparing r decays to one and 
three charged tracks and with or without a neutral pion. 
Uncertainties introduced by the particle identification for 
kaons and leptons are studied using control samples with 
high purities for the particles in question. The impact of 
the finite statistics of the simulated signal events is de- 
duced from the fit error of the efficiency determination. 

The uncertainty on the number of charged and neutral 
B mesons in the data set is determined as in [Hill and 
the branching fractions of the decays of the D* and the 
D are taken from [lit . 

Uncertainties introduced by the physics model which 



TABLE L Extracted yields for the four signal modes in the 
five relevant Am-spectra. 



mode 1 cos -dh 




Df 


Dt 


d;+ 


D*n+ [0.00|0.25] 


344 


273 


212 


152 


D*n+ [0.25|0.50] 


470 


238 


286 


123 


D*Tv+ [0.50|0.75] 


699 


170 


439 


83 


i3*7r+ [0.75|1.00] 


1027 


67 


668 


31 


Dn+ 




8414 




3361 



was used to simulate the MC data have been addressed 
by re-weighting the signal MC to an alternative decay 
model based on HQET The fit was repeated with 
efficiencies deduced from the reweighted signal MC and 
the deviations in the results are taken as systematic un- 
certainties. A possible influence of the background de- 
scription has been tested by varying the parametriza- 
tions. The backgrounds are alternatively described by a 
square root function, /(Am) = ^ Am — mo, where mg is 
the kinematic limit, multiplied by either polynomials or 
exponentials in Am. 

Table |TT] gives a summary of the various sources of 



TABLE II: Summary of systematic uncertainties of the deter- 
mination of the semileptonic branching fractions. 



Source 


AB{B - 


D**iv)/B{B - 


>D"iv)\%] 




D'l 


Df 


Dt 


DI+ 


tracking 


1.76 


1.39 


1.03 


1.14 


7r° efficiency 


0.06 


0.29 


3.25 


0.60 


particle identification 


2.61 


2.75 


3.11 


1.60 


MC statistics 


1.80 


5.61 


2.50 


3.32 


helicity correction 


0.65 


0.14 


0.17 


0.31 


number of B mesons 


2.68 


2.68 


2.68 


2.68 


B(D*+ D'\+) 


0.76 


0.19 


0.04 


0.10 




0.11 


0.45 


5.07 


0.93 


B{D° K-TV+) 


1.89 


0.42 


1.78 


2.03 


B{D+ K-TV + TT+) 


0.07 


2.67 


0.24 


0.54 


signal modeling 


2.11 


4.75 


3.21 


1.95 


bkg. parametrization 


1.93 


1.68 


3.20 


2.71 


total 


5.76 


9.03 


9.16 


6.17 
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systematic uncertainty and their impact on the results. 
Added in quadrature the total systematic uncertainties 
in the semileptonic branching fractions are 6-10%, de- 
pending on the D** type. 



In summary, we have measured the four branching 
fractions of B mesons decaying semileptonically into nar- 
row D** states. The D** decay rates are unknown, thus 
we can only determine the product branching fractions: 



B{B+ D\tvi) X B{D\ D*+TT-) 

B{B° D-e+i^e) X ^(L'f ^ 
B{B° D;-£+ue) x B{D;- D^*^"n-) 



= (2.97 ± 0.17,tat ± 0.17syst) X 10-3, 

= (2.29 ± 0.23stat ± 0.21syst) x W^, 

= (2.78 ± 0.24,tat ± 0.25,y,t) x 10~\ 

= (1.77 ± 0.26stat ± 0.11.y.t) X 10-3. 



We observe all modes with significance greater than 5f7, 
among them evidence of the Dj" contribution to the de- 
cay B D*Tr£i'. For modes already observed we find 
results in agreement with previous measurements, but 
achieve better precisions 0,1^, [l3| . 

For the decays of the D** we measure the branching 
ratio Bu/jji*) = 0.62 ± O.OSstat ± 0.02syst- This ratio is in 
agreement with theoretical predictions and previous 
measurements fill but has a smaller uncertainty by a 
factor of about four. 

For the Di we determine the polarization parameter to 
he Ad-^ = 3.8±0.6stat±0.8syst- It is the first measurement 
of the Di polarization, within the uncertainties consistent 
with unpolarized Di decaying purely via D wave, which 
gives the prediction =3, but violates HQS [J|. 
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